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CHAPTER  I 
INTRODUCTION 

In  the  Natural  Gas  Industry  production  rate  decisions  are 
dependent  upon  many  factors.  Ideally,  all  of  these  should  be 
considered  when  determining  final  production  rates  for  individual 
we 1 1 s . 

For  example,  it  is  desirable  to  attempt  to  produce  wells 
rateably.  That  is,  to  produce  more  from  the  underproduced  wells 
and  less  from  the  overproduced  wells.  The  final  goal  of  this 
method  of  determining  production  rates  is  to  have  all  wells 
depleted  at  the  same  point  in  the  future. 

At  the  same  time,  it  is  desirable  to  minimize  take-or-pay. 
Take-or-pay  is  the  penalty  resulting  from  the  underproduction  of 
a  well.  The  amount  of  take-or-pay  depends  upon  contractual 
agreements  between  the  producer  and  the  owner  of  the  well.  In 
addition,  some  amount  of  take-or-pay  can  usually  be  recovered  by 
overproduction  at  some  future  time.  This  again  depends  upon 
contractual  agreements.  Since  take-or-pay  is,  in  effect,  paying 
for  gas  not  received,  its  minimization  is  very  important. 

Another  consideration  to  be  made  when  determining  production 
rates  for  sources  is  system  capacity.  It  is  necessary  to 
consider  the  range  of  flow  rates  that  each  part  of  the 
transportation  system  is  physically  capable  of  handling. 

Minimizing  take-or-pay  may  require  a  maximum  production  rate  from 
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every  source  in  one  field,  while  requiring  a  minimum  production 
rate  from  every  source  in  another  field.  It  is  possible  that 
these  rates  will  not  be  within  the  acceptable  range  of  operating 
rates   for  the   available   transportation   system. 

Finally,  because  the  Natural  Gas  Industry  is  highly 
competitive,  it  is  important  to  minimize  cost.  Minimization  is 
necessary  if  a  competitive  price  for  natural  gas  is  to  be 
maintained. 

Each  of  the  factors  discussed  thus  far,  when  considered 
individually,  will  give  different  answers.  Ultimately,  a 
weighted  solution  is  desired.  The  first  step  in  the  development 
of  such  a  solution  involves  the  development  of  methods  for 
handling  each  of  the   factors   individually. 

To  further  improve  this  weighted  method  of  flow  rate 
determination,  it  is  desirable  to  consider  that  the  flow  rate 
from  each  source  may  be  changed  monthly,  and  that  a  desirable 
solution  over  several  time  steps,  or  months,  is  necessary.  In 
other  words  this  solution  should  be  applicable  not  only  to  one 
month  studies,  but  to  studies  covering  longer  time  periods  as 
well.  This  type  of  temporal  analysis  can  be  used  for  planning 
future    production    programs. 

As  a  first  step  in  the  solution  of  this  very  large  problem, 
this  paper  presents  a  method  of  determining  flow  rates  for  a 
large  number  of  natural  gas  sources  over  as  many  as  sixty  time 
steps.  The  program  developed  determines  the  flow  rates  for  a 
minimum  final    cost   of   natural    gas   taking  into  consideration: 

l)the   time  value   of   money 


2)  the  changing  prices  of  gas  from  existing  wells 

3) the  decrease  in  reserves  and  maximum  flow  rates  as  a 

result  of  gas  taken  from  a  well 
4)  the  need  for  each  well  to  have  the  ability  to  have  a 

unique  minimum  flow  rate 
5)the  need  for  the  user  of  the  program  to  have  the  ability 

to  change  the  desired  rate  for  each  time  period. 

Chapter  II  of  this  thesis  is  devoted  to  a  discussion  of  the 
literature  review  that  was  undertaken.   Nearly  all  of  this  review 
involved  the  use  of  the  computerized  search  service  available 
through  Farrell  Library.   Four  of  the  data  bases  which  were 
searched  are  discussed. 

In  Chapter  III,  a  discussion  of  the  logic  used  to  develop 
the  final  program  can  be  found.  This  discussion  carries  the 
reader  through  the  entire  program,  describing  each  of  the 
inportant  points. 

Finally,  Chapter  IV  presents  the  final  results  of  this 
research  project.  Several  special  data  sets  are  introduced,  run 
through  the  program,  and  the  results  discussed. 


CHAPTER  II 

LITERATURE  REVIEW 

At  the  onset  of  this  project,  the  author  of  this  paper  had 
no  feel  for  the  amount  of  work  that  had  already  been  done  in  this 
area.  To  insure  that  a  solution  had  not  already  been  published, 
an  extensive  literature  search  was  undertaken.  The  majority  of 
the  search  was  done  using  the  Computerized  Search  Service 
available  through  Farrell  Library.  A  discussion  of  the  results 
of  these  searches  follows. 

COMPENDEX  DATABASE 

Initially,  a  search  was  made  of  the  COMPENDEX  database. 
This  database  is  the  machine  readable  version  of  the  Engineering 
Index.  It  covers  approximately  3500  journals,  engineering  society 
and  organization  publications,  conference  proceedings  and 
government  reports.   Using  key  words  of 

1)  Extremum  of  Function  or  Optimize  or  Maximize  or 

Minimize  and 

2)  Natural  Gas  or  Well  and 

3)  Production 

resulted  in  51  abstracts  that  appeared  at  first  glance  to  be 
related  to  the  problem  under  study.  Of  these  abstracts,  only  one 
proved  to  be  of  any  value.  This  article,  entitled  "Computerized 
System  to  Optimize  Daily  Oil  and  Gas  Production  in  Kuwait"  by 
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Cain  and  Shehata  appeared  by  title  to  be  precisely  the  type  of 
article  desired.  Acquisition  of  the  article  proved  discouraging 
in  that  it  covered  the  optimization  of  scheduling  regular 
maintenance  with  production  schedules.  Although  not  the  same 
problem  as  the  one  under  study,  the  solution  arrived  at  by  Cain 
and  Shehata  was  one  requiring  the  use  of  linear  programming 
techniques.  Theirs  was  a  much  smaller  scale  problem  than  the  one 
considered  here.  They  state  that  "The  full-size  problem  normally 
consists  of  400  to  500  constraints  and  between  1,300  and  1,500 
variables."  Further,  they  indicate  that  the  solution  process  may 
be  rather  time  consuming  "In  operation,  1  to  1  1/2  hours  is 
allowed  for  running  the  model  and  dispatching  the  results  to  the 
field".  It  is  impossible  to  know  how  much  of  the  time  alloted 
was  necessary  for  the  problem  solution  and  how  much  was  for 
dispatching  the  results.  However,  even  if  only  30  minutes  of  the 
time  was  required  for  problem  solution,  it  becomes  evident  that 
a  problem  of  the  magnitude  being  considered  in  this  thesis  could 
not  be  solved  in  the  same  manner  without  requiring  an  excessive 
amount  of  computer  time. 

As  is  often  the  case,  the  paper  by  Cain  and  Shehata  contained 
a  reference  in  its  Bibliography  to  another  publication  that  was 
more  helpful  than  the  original.  This  publication  was  entitled 
"Mathematical  Models  to  Help  Manage  the  Oil  Resources  of  Kuwait" 
by  Ali,  Beasley,  Batchelor,  and  Beale.  This  article  proved  to  be 
more  along  the  lines  of  the  problem  being  considered  here. 
Although  the  problem  was  not  exactly  like  that  being  studied, 
there  were  enough  similarities  to  confirm  that  the  method  of 
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solution  being  used  was  an  accepted  method.  The  Ali,  Beasley, 
Batchelor,  and  Beale  article  states:  "  The  first  joint 
development  was  a  model  to  specify  a  schedule  of  wells  to  use  in 
meeting  a  given  daily  total  production  rate.  The  schedule  has  to 
operate  within  the  capacity  limits  of  the  plant,  produce  oil  of 
the  right  quality,  satisfy  demands  on  the  various  gas  streams  and 
optimize  liquid  gas  production  while  attempting  to  meet  reservoir 
engineering  requirements:  This  is  a  natural  application  of 
linear  programming,...".  Even  though  this  description  sounds 
very  much  like  the  problem  under  study,  there  was  one  very 
important  difference,  that  being  again,  the  magnitude  of  the 
problem.  The  Ali,  Beasley,  Batchelor,  and  Beale  article  stated, 
near  it's  conclusion  "  The  latest  problem  solved  has  2336 
constraints  and  4392  linear  programming  variables.  So. ..each  new 
problem  still  requires  several  hours  of  CPU  time  on  the  Univac 
1108  Computer."  If  4392  variables  required  several  hours  of  CPU 
time,  the  time  required  by  300,000  variables  was  sure  to  be 
prohibative. 

At  this  point  in  the  search,  two  articles  had  been  found 
that  described  problems  similar  to  that  being  studied,  both  used 
typical  linear  programming  techniques,  and  both  indicated  that  a 
considerable  amount  of  computer  time  was  required  for  a  small 
number  of  variables. 

INSPEC  DATABASE 

Another  database  known  as  the  INSPEC  was  searched.  This 
database  covers  the  fields  of  physics,   e  lectrotechno logy , 


computers,  and  control.  This  database  includes:  journal  papers, 
conference  proceedings,  technical  reports,  books,  patents,  and 
university  thesis.    Using  key  words  of 

1)  Gas  and 

2)  Natural  Gas  or  Well  and 

3)  Automatic  Control  or  Computer  Control  or  Computer 

Optimization  or  Computer  Production 

resulted  in  40  abstracts.  Many  of  these  dealt  with  topics  much 
different  than  the  one  being  studied.  Some  of  the  areas  brought 
out  by  this  search  included: 

1)  Geothermal  Wells 

2)  Microprocessor  based  controls 

3)  Automatic  Pipeline  Control  and 

4)  Chemical  Controls 

None  of  these  were  of  any  value.  However,  several  did  mention 
optimization  problems  that  required  the  use  of  linear  programming 
techniques. 

DOE  ENERGY  DATABASE 
Finally,  the  DOE  ENERGY  database  was  searched,  and  abstracts 
were  received.  This  database  is  one  of  the  worlds  largest 
sources  of  literature  references  on  all  aspects  of  energy  and 
related  topics.  It  provides  coverage  of  journal  articles,  report 
literature,  conference  papers,  books,  patents,  dissertations  and 
translations.   Using  the  key  words  of 

1)  Production  and 

2)  Natural  Gas  or  Well  and 

3)  Computer  Optimization 
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resulted  in  34  abstracts.  Of  these  34,  two  were  those  already 
discussed,  namely  the  Cain  and  Shehata  article,  and  the  Ali, 
Beasley,  Batchelor,  and  Beale  article.  Additionally,  an  article 
entitled  "Seeking  Optimum-Profit  Production  Decisions"  by  Cavaroc 
and  Sylla  was  discovered.  This  article  discussed  a  solution  to  a 
problem  very  similar  to  the  one  being  studied.  It  summarized  the 
problem  as  one  which  requires  the  "determination  of  production 
rates  that  will  yeild  the  maximum  profit  from  this  reservoir, 
subject  to  the  unit  allowable  and  production  characteristic 
limitations."  Although  never  stated,  the  problem  discussed  by 
Cavaroc  and  Sylla  was  probably  much  smaller  in  scope  than  the  one 
to  which  this  thesis  addresses  itself.  Several  places,  the  use 
of  a  matrix  solution  was  mentioned.  The  authors  stated:  "Another 
advantage  of  this  approach  is  the  simplicity  of  the  mathematics 
involved,  which  consists  merely  of  a  logical  and  systematic 
algebraic  expression  of  the  problem,  supplemented  with  a  limited 
knowledge  of  matrix  algebra."  A  matrix  solution,  although  very 
common,  is  unreasonable  to  use  for  a  problem  of  the  size  being 
considered  here.  It's  use  by  Cavaroc  and  Sylla  only  strengthens 
the  belief  that  a  solution  to  a  problem  of  this  magnitude  had  not 
been    published. 

APALIT  DATABASE 
A  short  search  through  APALIT  database  was  conducted.  This 
data  base  covers  all  aspects  of  information  dealing  with  the 
American  Petroleum  Institute.  It  primarily  covers  journal 
articles,  conference  proceedings,  books,  and  reports  dealing  with 
petroleum    and    petroleum    products    from    the    wellhead    through    the 
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pipeline.  This  search  revealed  so  little  that  a  copy  of 
abstracts  was  not  even  requested.  Although  this  search  produced 
nothing  of  any  value  the  fact  that  it  was  conducted,  and  that  the 
American  Petroleum  Institute  database  gave  no  worthwhile 
references    is   worth  mentioning. 

In  addition  to  the  computerized  literature  searches,  many 
text  and   other   reference   books  were  used   to  complete    this    thesis. 

The  literature  search  conducted  as  part  of  the  completion  of 
this  project  has  been  informative.  It  has  indicated  that  this 
specific  problem,  in  this  magnitude  has  not  previously  been 
considered.  Algorithms  exist  which  handle  much  smaller  sets  of 
data.  However,  for  this  problem,  the  increase  in  problem  size 
brought  on  the  realization  that  previously  used  techniques  of 
matrix   solutions  were  not  applicable. 


CHAPTER  III 

METHODS  AND  PROCEDURES 

A  set  of  equations  were  developed  which  describe  the 
relationship  between  the  amount  of  gas  taken  from  each  source  in 
each  time  step  and  the  overall  gas  cost.  These  equations  were 
supplemented  by  a  large  number  of  constraints  that  limited  the 
solution  to  the  problem.  Since  the  problem  was  a  linear  one, 
linear  programming  techniques  were  investigated,  and  were 
partially  used  in  the  development  of  the  final  computer  program. 
The  steps  followed  from  initial  equation  development  to  final 
program  results  will  be  detailed  in  this  chapter. 

Theoretical 
Cost  Equation 

Linear  problems  that  need  to  be  maximized  or  minimized  must 
first  be  formulated  into  an  equation.  This  equation  must  relate 
all  of  the  unknown  variables  to  the  parameter  being  maximized  or 
minimized.  In  this  case,  an  equation  needed  to  be  developed  that 
related  the  total  cost  of  the  gas  produced  over  "m"  time  periods 
to  the  rate  of  gas  production  from  "n"  sources  in  each  of  those 
time  periods. 

Several  considerations  were  made  in  the  development  of  the 
cost  equation  to  make  it  more  realistic.  The  equation  in  its 
simplest  form  is  shown  below. 
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Cost=R1/1(Clfl)+Rlf2(C1/2)+R1/3(C1/3)+----+ 
Rl,n(Cl,n)+R2,l(C2,l)+R2/2(C2,2)+ 

R2,3(C2,3)+ +R2,n(C2.n)+ + 

^,l<cm,l>+V2<cin,2>+Rm,3<C!n,3>+----  + 

Rm,n  ^Cm,n' 

Where  R-rate  of  flow  from  any  given  source 
C-cost  of  the  gas 

m-total  number  of  time  periods  being  considered 
n-total  number  of  sources  being  considered 

The  first  subscript  on  each  of  the  terms  refers  to  the  time 

period  being  considered  (up  to  60),  while  the  second  subscript 

refers  to  the  source  being  considered  (up  to  5000). 

In  order  to  make  the  individual  cost  terms  more  realistic, 
several  factors  were  introduced.  Each  gas  source  was  assigned  an 
initial  gas  cost  for  time  period  one.  This  cost  was  then 
adjusted  for  each  subsequent  time  period,  taking  into 
consideration  the  time  value  of  money,  and  the  actual  escalation 
of  natural  gas  prices. 

The  magnitude  of  the  escalation  factor  was  dependent  upon 

the  Natural  Gas  Policy  Act  (NGPA)  code  assigned  to  the 

source,  and  in  some  cases  upon  the  month  of  the  year.  Each  of 

these  escalated  prices  were  then  brought  back  to  present  worth 

using  the  equation: 

P  =  F  [1/  (l  +  i)n] 

Where  P-  present  worth  of  the  cost  term 
F-  future  worth  of  the  cost  term 
i-  inflation  rate  per  time  period 
n-  number  of  time  periods 

The  individual  cost  terms  in  the  overall  cost  equation  can  be 

mathematically  described  by: 
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Cm  ri=C1  _  (escalation  factor) m_1 
(inflation  factor) m~1 


Constraints 

Earlier,  introduction  of  constraints  used  to  limit  the 
problem  were  mentioned.  One  of  these  constraints  was  a  specified 
flow  rate  to  be  achieved  in  each  time  period.  The  sum  of  the 
rates  from  each  source  in  a  given  time  period  must  be  equal  to 
this  specified  rate.  The  user  of  the  program  must  specify  these 
desired  rates. 

Once  specified,  these  rates  were  labeled  R(l),  R(2), 
R(3)...R(m).   In  equation  form 

R(l)=R1/1+R1/2+R1/3+....+Rlfn 

R(2)=R2/1  +  R2/2  +  R2,3  +  ....+R2fI1 

R<m>=Rm,l+Rm,2+Rm,3  +  —-+Rm,n 
These    desired    rate    equations    were    used    to    eliminate    the    highest 

priced   gas    source    from   the   cost   equation.      The   highest   priced 

source  was  the  one  chosen  to  be  eliminated  for  several   reasons. 

First,    since    the   cost  was    to   be   minimized,    the   highest   priced 

source  was   the    least    likely   to   be      assigned   a   gas   flow   above    its 

minimum  required  flow.      Second,   eliminating  the  highest  priced 

source  resulted  in  all    negative   coefficients    in  the   overall    cost 

equation.       The    need    for    negative    coefficients    will    become    clear 

later. 

The  desired  rate  equations  were  rewritten  to  isolate  the 

highest  priced  source.   If  the  subscript  "k"  is  used  to  indicate 


13 


the  highest  priced  source  then 

Rl,k=R(1)"Rl,l"Rl,2~Rl,3--  •"Rl,Jc-l~Rl,k  +  l*  '  '"Rl,n 

R2 , k=R ( 2  >  ~R2 , 1"R2 , 2"R2 , 3  * ' ' "R2 ,k-l"R2 , k+ 1  *  * • "R2 , n 

Rm>k=R(m)  -Rm,i-Rm,2"Rm,3*  •  •~Rm,k-l'_Rra,k+l'  *  '"^n 

Substituting  these  into  the  cost  equation  results  in  the  equation 
below. 

Cost=R1/1(Clfl)+R1/2(Clf2)+R1>3(C1/3) 

♦•••+Rl,k-l<Cl,k-l)+lR(1)~Rl,l"Rl,2 
"Rl,3"*~Rl,k-l~Rl,k+l* ' -"Rl,n]cl,k 

+Rl,k+l(cl,k+l)+,,-+Rl,n(cl,n)+- •• 

R2,1(C2,1)+R2,2(C2,2)+R2,3(C2,3) 

+---+R2,k-l(C2,k-l)+[R(2)-R2,l-R2/2 
-R2/3.  ••-R2,k-l"R2,k+l","R2,n]C2,k 
+R2,k+l(C2,k+l)+'  ••+R2,n(C2,n)  +"  • 

Rm,l<Cm,l>+Rm,2(Cm,2>+V3<Cm,3> 

+---Rm,k-l<Cm,k-l)+tR(m)-Rm/l-Rm,2- 

~Rm,3* ' •"Rm/k-l"Rm,k+l- ' •~Rm,nJcm,k 

+Rm,k+l(cm,k+l) +* ■ * +Rm,n (cm,n) 
Regrouping  terms  gives  the  equation  below. 

Cost=R(l)C1  k+R(2)C2/k+R(3)C3fk+. . .+R(m)Cm^k+ 

Rl,l(Cl,l-Cl,k>+Rl,2<Cl,2-Cl,k>+Rl,3<Cl,3-Cl,k>+'--+ 
Rl,k-l(cl,k-l"cl,k)+Rl,k+l(cl,k+l"cl,k) + • ■ ■ +R1 ,n (cl ,n~ 
Cl,k)+R2,l(C2,l"C2,k)+R2,2(C2,2"C2,k) +R2,3 (C2,3" 
C2,k)  +  "-+R2,k-l(c2,k^l~c2,k)  +R2,k+1  (c2  ,k+l~c2  ,k'  +-  •  -  + 
R2,n(c2,n-c2,k'+-"+Rm,l(Cm,l-Cm,k)+Rm,2(Cm,2-Cm,k)  + 
Rm,3 (cm,3-cm,k) +- • -+Rm,k-1 (Cm,k-1-Cm,k) + 
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Rm,k+l(Cm,k+l~Cm,k)+* ' • +Rm,n (cm,n"Cm,k) 
This  substitution  has  resulted  in  the  cost  equation  having 

several  interesting  and  helpful  characteristics.  First,  the  first 

"m"  terms  represent  the  maximum  possible  cost  of  the  gas.   In 

other  words,  these  represent  the  cost  of  the  gas  if,  in  each  time 

step,  the  gas  was  all  taken  from  the  highest  priced  source.   This 

is  a  worst  case  situation  that  can  only  be  improved.   Second, 

each  of  the  rate  coefficients  has   taken  on  a  negative  value. 

Each  of  these  coefficients  is  the  difference  between  the  cost  of 

gas  for  a  given  source  in   a  given  time  period  and  the  cost  of 

gas  from  the  highest  priced  source  in  the  same  time  period. 

Clearly,  upon  reflection,  to  minimize  the  overall  cost,  or  to 

decrease  the  large  positive  number  at  the  beginning  of  the  cost 

equation,  it  will  be  necessary  to  take  as  much  gas  as  possible 

from  the  source  with  the  next  largest  negative  coefficient.   This 

process  should  be  repeated  until  all  flow  requirements  are  met. 

The  introduction  of  a  limited  production  range  for  each 
source  is  another  constraint  used  to  contain  the  problem.  The 
user  must  specify  an  initial  minimum  and  maximum  allowable 
production  rate.  The  minimum  rate  remains  constant  throughout 
the  study  while  the  maximum  rate  decreases  as  the  supply  of 
natural  gas  available  decreases. 
Solution  Format 

Problems  similar  to  the  problem  under  study  are  usually 
solved  through  the  use  of  linear  programming  techniques. 
However,  such  an  approach  gives  rise  to  a  matrix  solution.  In 
addition,  this  appproach  requires  the  introdution  of  slack 
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variables  into  the  problem.  Slack  variables  are  variables  used 
to  remove  inequalities  from  restraint  equations.  For  example, 
the  equation 

would  become,  with  the  introduction  of  slack  variables,  the 
following  two  equations 

R1,1*R2,1+R3,1+R4,1+S1"R(1) 

0<S1<R(1) 

One  such  variable  is  required  for  each  constraint  that  involves 
an  inequality.  Since  the  rate  of  gas  flow  from  each  source  is 
constrained  by  a  maximum  and  a  minimum  desired  flow  rate,  these 
constraints  are  all  inequalities.  As  a  result,  one  slack 
variable  would  have  to  be  introduced  for  each  unknown  rate. 
Since  the  problem  at  hand  could  involve  up  to  300,000  unknowns 
(5000  sources  *  60  time  steps)  before  any  slack  variables  are 
introduced,  a  matrix  solution  seemed  impractical,  hence  a 
solution  requiring  less  active  storage  space  on  the  computer  was 
desired. 
Computer   Programs 

The  solution  used  required  the  development  of  a  series  of 
two  computer  programs  capable  of  handling  a  problem  containing  up 
to  300,000  unknown  variables  without  using  an  excessive  amount  of 
active  computer  memory.  Copies  of  these  programs  can  be  found  in 
Appendix  I.  The  final  programs  developed  use  two  data  files  to 
handle  all  of  the  necessary  information.  The  first  data  file 
contains  one  card  for  each  desired  rate  (up  to  300,000  cards) 
while  the  second  data  file  contains  one  card  for  each  source   (  up 
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to  5000  cards) . 

The  desired  rate  cards  contain  four  pieces  of  data:  namely 

1)  The  source  number 

2)  The  time  step  being  considered 

3)  The  coefficient  for  the  rate  term  in  the 

cost  equation 

4)  The  desired  rate  (the  unknown) 

These  cards  are  generated  by  the  program  rather  than  being 
supplied  by  the  user. 

The  source  cards  contain  nine  pieces  of  data:  namely 

1)  The  state  code 

2)  The  field  code 

3)  The  source  name  (well  name) 

4)  A  flag  to  indicate  whether  the  well  is  deep 

or  shallow 

5)  The  allowable  production  rate  for  the  source 

6)  The  minimum  production  rate 

7)  The  maximum  production  rate 

8)  The  price  of  the  gas  in  time  period  one 

9)  A  flag  to  indicate  whether  the  NGPA  code 

is  102  or  108,  or  not  102  or  108 

The  information  required  for  these  cards  is  supplied  by  the  user. 
The  desired  rate  cards  are  created  and  sorted.  Following 
this,  they  are  accessed  only  once.  The  source  cards,  on  the 
other  hand,  are  created  and  then  referenced  up  to  60  times  each 
in  an  unpredictable  order.  As  a  result,  the  source  card  file  is 
a  direct  access  input/output  file  (Direct  Access  I/O).  The 
desired  rate  file  is  a  sequential  file. 
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Each  of  the  pieces  of  data  on  these  two  data  cards  deserves 
particular  discussion. 

The  first  two  items  on  the  desired  rate  cards  fully  describe 
the  case  being  determined.  For  example,  one  card  may  refer  to 
the  67th  source  in  the  23rd  time  period,  while  another  may  refer 
to  the  68th  source  in  the  23rd  time  period.  The  third  item  is 
the  actual  coefficient  as  it  would  be  determined  for  the  cost 
equation.  This  would  be  the  coefficient  that  would  appear  in  the 
final  form  of  the  cost  equation.  The  last  item  is  the  desired 
output  of  the  program,  the  answer. 

The  first  three  items  on  the  source  cards  fully  identify  the 
gas  source.  The  cards  are  in  a  particular  order,  and  the  source 
number  used  on  the  desired  rate  cards  corresponds  to  the  position 
in  the  data  file  of  the  source  card.  For  example,  the  first 
source  card  would  have  a  source  number  of  one,  the  third  a  source 
number  of  three  etc.  As  a  result,  it  is  not  necessary  to  have 
the  source  number  as  one  of  the  identifying  pieces  of  data  on  the 
source  card.  The  fourth  item  on  the  source  card  is  a  flag 
indicating  the  relative  depth  of  the  source.  This  piece  of 
information  is  used  to  calculate  the  initial  reserves  for  a  given 
source.  A  shallow  well  (depth  indicator  =  1)  is  assumed  to  have 
initial  reserves  equal  to  20  years  production  at  the  monthly 
allowable  rate  or  240*A1  lowable. 

A  deep  well  (depth  indicator  =  0)  is  assumed  to  have  initial 
reserves  equal  to  10  years  production  at  the  monthly  allowable 
rate  or  120*A1 lowable. 

The  fifth  item  on  the  source  card  is  the  allowable  monthly 
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production  rate  for  the  well.   This   value  is  used  along  with  the 
depth  indicator  to  calculate  the  initial  reserves  for  the  source. 

The  sixth  item  is  the  minimum  production  rate  for  the  well. 
This  value  corresponds  to  the  smallest  amount  of  gas  that  can  be 
taken  from  a  source  in  any  given  time  period  without  violating 
contractual  agreements  between  the  producer  and  the  company 
buying  the  gas. 

The  seventh  item  is  the  maximum  production  rate.  This 
maximum  rate  decreases  as  gas  is  taken  from  a  source.  This 
decrease  is  calculated  by  finding  the  percentage  decrease  in  the 
reserves,  and  applying  that  same  percentage  decrease  to  the 
maximum  rate. 

The  eighth  item  is  the  initial  price  of  the  gas.  The 
calculation  of  the  coefficient  for  the  cost  equation  depends 
upon  this  value. 

The  last  item  is  the  NGPA  code  flag.  All  wells  are 
escalated  in  price  .5%  per  month.  In  addition,  wells  with  an 
NGPA  code  of  102  or  108  (NGPA  code  flag=l)  are  escalated  in  price 
an  additional  4%  each  December.  These  considerations  go  into  the 
calculation  of  the  cost  equation  coefficients. 

Two  flow  charts  showing  the  basic  steps  followed  by  the 
computer  programs  in  solving  the  problem  are  shown  the  Appendix 
II.  As  shown  in  the  first  flow  chart,  the  user  must  supply 
several  pieces  of  data.   These  include: 

1)  The  number  of  sources  being  considered 

2)  The  number  of  time  steps  being  considered 

3)  The  beginning  month  (1-12  January-December) 
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4)  The  monthly  inflation  factor  to  be  used 

5)  The  desired  rates  for  each  of  the  time 

steps 

6)  The  information  for  the  source  cards. 

Once  items  one  through  five  above  have  been  accepted  by  the 
program,  the  first  source  card  is  read  in.  The  allowable  rate 
and  the  depth  indicator  from  that  source  card  are  used  to 
calculate  the  initial  reserves  for  the  source.  The  source  is 
then  tested  to  see  if  it  is  the  highest  priced  source  encountered 
thus  far.  If  it  is,  the  source  number  (  the  number  of  source 
cards  thus  far  read  in)  is  saved,  along  with  the  price  and  the 
NGPA  code  associated  with  that  source.  As  a  final  step,  before 
the  source  card  is  created,  the  minimum  rate  for  this  source  is 
added  to  a  running  total  of  all  of  the  minimum's  for  all  of  the 
sources  being  studied.  Finally,  a  card  in  the  source  data  file 
is  created  using  the  information  from  this  source.  This  cycle  is 
repeated  once  for  each  source  in  the  study. 

When  completed,  all  of  the  data  cards  in  the  source  file  have 
been  created,  the  highest  priced  source  has  been  identified,  and 
the  algebraic  sum  of  all  of  the  required  minimum  flow  rates  has 
been  calculated.  The  source  number,  initial  price  and  NGPA  code 
for  the  highest  priced  source  have  also  been  identified  and 
saved. 

Once  the  source  file  has  been  created,  the  program  begins 
the  process  of  creating  the  desired  rate  file.  The  source  file 
is  referenced  one  element  at  a  time,  beginning  with  the  first 
element  and  proceeding  sequentially  through  the  file.   A  month 
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flag  is  set  to  the  initial  month  (specified  by  the  user  at  the 
beginning  of  the  program)  ,  and  the  Price  and  NGPA  code  are  read 
off  of  the  source  card.  A  check  is  made  to  insure  that  the 
source  being  considered  is  not  the  highest  priced  source.  If 
such  is  the  case,  the  program  skips  out  of  the  loop,  and  selects 
the  next  source  card.  Finally,  the  coefficient  for  the  source 
being  considered  is  calculated  for  this  time  period  using  the 
initial  price,  the  escalation  factor  (based  on  NGPA  code),  the 
inflation  factor  and  current  time  period  for  the  source  being 
considered  as  well  as  for  the  highest  priced  source.  This 
coefficient  is  then  used  to  create  a  data  file  element. 
Additionally,  an  initial  rate  equal  to  the  minimum  required  rate 
is  written  to  the  data  card.  The  desired  rate  card  contains,  at 
this  point,  the  source  number,  the  time  step  number,  the 
calculated  coefficient,  and  the  desired  rate  (initially  set  to 
the  minimum  required  rate). 

The  steps  described  above  are  repeated  once  for  each  time 
step  to  be  considered  without  changing  the  source.  Before  each 
pass  through  the  loop,  the  time  step  number  and  month  are 
incremented  by  one.  When  the  first  set  of  these  passes  is 
completed,  the  desired  rate  file  contains  a  number  of  data  cards 
equal  to  the  number  of  time  steps  being  considered.  Each  data 
card  has  the  same  source  number  and  the  same  desired  rate,  but 
different  time  step  numbers  and  different  coefficients. 

The  process,  beginning  with  the  reading  of  a  source  card,  is 
repeated  until  all  source  cards  have  been  read.  This,  when 
completed  will  have  resulted  in  the  creation  of  the  entire 
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desired  rate  file.  At  this  point,  the  file  contains  data  cards 
arranged  by  source  number  and  time  step,  with  all  desired  rates 
set  at  the  sources  minimum  acceptable  flow  rates. 

As  discussed  earlier,  all  of  the  coefficients  should  be 
negative  and  will  need  to  be  dealt  with  in  order  of  decreasing 
absolute  value.  The  program  deals  with  this  problem  by  using  a 
Sort  Processor  to  sort  the  data  cards  of  the  desired  rate  file. 
This  new  sorted  desired  rate  file  is  organized  in  such  a  manner 
that  its  elements  are  in  order  of  cost  equation  coefficients,  but 
can  be  linked  to  a  particular  source  and  time  step  by  the  first 
two  pieces  of  data  on  each  card. 

Before  the  second  program  deals  with  the  sorted  desired  rate 
file,  the  first  program  calculates  the  large  positive  number  at 
the  beginning  of  the  cost  equation  and  saves  that  information. 
This  number  represents  the  cost  of  the  gas  if  all  of  the  gas  for 
each  time  step  is  taken  from  the  most  expensive  source. 

Once  all  of  the  above  steps  have  been  completed,  the  second 
program  begins  execution,  and  the  first 

desired  rate  card  can  be  dealt  with.  This  card  should  represent 
the  source  and  time  step  with  the  largest  negative  coefficient, 
or  in  other  words,  the  least  expensive  gas  available.  This  card 
will  contain  information  relating  it  to  a  particular  source.  The 
source  card  for  that  source  must  also  be  used  in  the 
determination  of  a  final  desired  rate. 

Several  things  must  be  considered  when  determining  the  final 

desired  rate  for  a  source.   Namely: 

1)  Is  the  overall  desired  rate  for  this  time  step 
satisfied? • 
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2)  If  the  overall  desired  rate  is  not  satisfied,  is  it 

less  than  the  maximum  allowable  rate  for  this 
source? 

3)  Are  the  reserves  for  this  source  depleated? 

4)  If  the  reserves  are   not  depleated,  are  they  less 

than  the  maximum  allowable  rate  for  this  source? 

Taking  into  consideration  all  of  these  things,  the  final  desired 
rate  is  determined. 

Once  a  final  desired  rate  has  been  determined,  several  of 
the  parameters  of  the  problem  must  be  adjusted.   These  include: 

1)  The  desired  rate  for  the  time  step 

2)  The  reserves  of  the  source 

3)  The  maximum  allowable  rate  for  the  source 

The  last  two  of  these  changes  replace  previous  values,  and  are 
written  back  to  the  source  card  file  with  the  other  source 
information. 

As  a  final  set  of  steps  in  this  portion  of  the  program,  the 
final  results  are  printed  out,  and  the  total  cost  term  is 
adjusted  by  the  product  of  the  coefficient  of  this  source  and  its 
final  desired  rate.  This  product  will  be  a  negative  number  and 
will  decrease  the  large  positive  number  at  the  beginning  of  the 
cost  equation. 

The  entire  process  is  repeated  until  all  of  the  coefficient 
cards  have  been  examined  in  their  sorted  order.  At  this  point  in 
the  solution  process,  a  final  desired  rate  has  been  assigned  to 
all  of  the  sources  for  each  of  the  time  steps  with  the  exception 
of  the  highest  priced  source.  This  source  is  assigned  a  rate 
above  its  minimum  rate  only  if,  after  all  of  the  other  sources 
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have   been  assigned  final    rates,    the  desired  rate  for      any  time 
period   is   still    not   satisfied. 

Finally,    the  program  prints  out  the  final   overall   cost. 


CHAPTER  IV 


RESULTS  AND  DISCUSSION 


This  thesis  has  presented,  to  this  point,  the  logic  used 
and  steps  followed  in  the  development  of  a  computer  program  to 
minimize  the  cost  of  natural  gas  from  a  large  number  of  sources 
over  a  long  time  frame.  The  program  was  written  to  handle  as 
many  as  5000  gas  sources  over  as  many  as  60  time  steps. 

To  show  that  the  original  objective  of  the  project  has  been 
met,  several  special  sets  of  data  were  used.  The  results  of 
these  special  cases  will  be  discussed  in  this  chapter.  All  of 
the  numbers  used  in  the  data  sets  are  realistic  numbers,  but  have 
been  identified  by  source  names  that  cannot  be  traced  to  actual 
gas  sources.  This  has  been  done  to  protect  against  the  release 
of  any  proprietary  information  which  may  have  been  used  in  the 
development  of  this  program. 

The  first  data  set  used  consisted  of  seventy-one  gas  sources 
and  used  the  full  60  time  step  range.  A  copy  of  this  data  set, 
labeled  DATA  SET  I  can  be  found  in  Appendix  III.  This  was  the 
largest  set  of  data  used,  and  it  was  intended  to  show  several 
desired  results.  First  of  all,  the  program  was  handling,  with 
this  data  set,  4260  unknowns.  This  is  a  number  comparable  with 
several  of  the  data  sets  discussed  in  Chapter  II  .  The  total 
time  reguired  by  the  program  to  execute  when  dealing  with  this 
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number  of  unknowns  was  four  minutes  and  thirteen  seconds.  This 
is  a  very  short  execution  time  when  compared  with  any  of  the 
programs  discussed  in  Chapter  II.  Undoubtedly  the  amount  of  time 
required  to  execute  the  program  with  the  largest  possible  data 
set  would  be  very  large,  but  this  program  appears  to  be  much  more 
efficient  than  those  already  in  existence. 

A  second  result  that  can  be  shown  from  this  large  data  set 
is  that  the  desired  rate  specified  in  the  data  cards  is  met.  An 
entry  in  Appendix  III  labeled  DESIRED  RATE  SUM  FOR  ALL  FIELDS  IN 
TIME  PERIOD  ONE  (DATA  SET  I)  shows  this.  On  this  summary  sheet 
are  all  of  the  assigned  rates  for  the  first  time  period.  One 
rate  is  given  for  each  source  in  the  study.  Their  sum  is  the 
desired  rate  specified  in  the  input  data  set. 

A  second  data  set  almost  exactly  the  same  as  the  first  was 
used  to  show  that  the  program  does  indeed  minimize  cost.  In  this 
data  set,  the  least  expensive  gas  source  was  given  a  minimum  and 
a  maximum  flow  rate  of  zero.  This  should  have  forced  the  program 
to  satisfy  the  required  rates  by  using  more  expensive  gas  than 
would  have  been  necessary  with  the  first  data  set.  As  a  result, 
the  final  cost  of  the  gas  should  have  risen  slightly  to  indicate 
that  the  program  was  responding  as  expected.  The  total  cost  of 
the  gas  when  utilizing  the  cheapest  source  was  $393,810,044.00, 
while  the  cost  of  the  gas  when  the  cheapest  source  was  restrained 
was  $394,187,640.00.  An  entry  in  Appendix  III  entitled  "Final 
Costs  For  60  *  71  Data  Sets"  shows  these  final  costs.  This 
slight  increase  in  cost  was  precisely  what  was  expected. 

The.  last  three  data  sets  were  smaller  than  the  first  two, 
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and  were  used  to  show  that  cost  of  the  gas  was  calculated 
properly,  that  the  program  followed  given  constraints  and  that 
the  desired  rates  were  satisfied. 

The  first  of  these  smaller  data  sets  restrained  all  but  the 
cheapest  of  twelve  sources.  As  a  result  only  the  least  expensive 
source  should  have  been  assigned  a  flow  rate.  A  copy  of  this 
data  set  (labeled  DATA  SET  II),  and  part  of  the  results  printed 
out  by  the  program  (labeled  RESULTS  OF  DATA  SET  II)  are  shown  in 
Appendix  III. 

Similarly,  in  the  second  of  the  samller  data  sets,  all  of 
the  sources  except  an  intermediate  priced  well  were  restrained. 
This  should  have  resulted  in  only  the  intermediate  priced  well 
receiving  an  assigned  flow  rate.  Again,  the  results  proved  to  be 
what  was  expected.  Part  of  the  results  printed  out  by  the 
program  (labeled  Results  of  Data  Set  III)  are  shown  in  Appendix 
III. 

Finally,  the  last  of  the  smaller  data  sets  restrained  all 
but  the  most  expensive  of  the  twelve  sources,  and  had  the 
expected  results.  Part  of  these  results  are  shown  in  Appendix 
III  (labeled  Results  of  Data  Set  IV). 

Several  results  from  the  last  three  data  sets  discussed 
deserve  particular  attention. 

The  final  costs  resulting  from  these  data  sets  should  be  in 
an  expected  order,  with  the  least  expensive  source  producing  the 
lowest  final  cost,  the  intermediate  source  producing  an 
intermediate  final  cost,  and  the  highest  priced  source  producing 
the  highest  final  cost.  Copies  of  the  computer  printout  showing 
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these  total  costs  can  be  found  in  Appendix  III  (labled  Final 

Costs  for  Data  Sets  II,  III,  and  IV).   The  actual  results  are 

shown  in  the  table  below. 

TOTAL  COST  TABLE  FOR  DATA  SETS  II, III,  AND  IV 

Description  of  Data  Total  Final  Cost 

Data  Set  II  $11,023.98 

Data  Set  III  $72,959.63 

Data  Set  IV  $105,023.67 

These  results  are  exactly  what  was  expected. 

Also  worthy  of  mention  is  the  match  between  the  actual 

assigned  rates  for  the  sources  and  the  desired  rates  specified  by 

the  input  data  set.   For  example,  Data  set  II  (Appendix  III) 

3 
shows  a  desired  rate  for  time  period  three  as  1.8  *  10  .   The 

output  for  Data  set  II  (Appendix  III)  shows  time  step  three  with 

an  assigned  rate  of  1800  units. 

Finally,  the  order  of  the  assigned  rates  on  the  printout  is 
of  particular  significance.  Since,  in  this  particular  data  set, 
(Data  set  II)  the  starting  month  was  specified  as  month  4,  and 
since  it  is  being  compared  with  the  highest  priced  source  for 
coefficient  calculations,  the  ninth  month  or  December  has  the 
largest  coefficient.  This  results  because  the  highest  priced 
source  is  one  that  is  escalated  an  extra  4%  in  December. 

As  a  last  check,  sample  calculations  for  the  determination 
of  the  maximum  flow  rate  reduction,  and  for  the  calculation  of  a 
coefficient  are  included  in  Appendix  IV.  The  computer  printout 
showing  the  actual  maximum  flow  rate  reduction  (labeled  Data 
Supporting  Maximum  Flow  Rate  Reduction  Calculations)  can  be  found 
in  Appendix  IV.   Similarly,  the  actual  coefficient  given  by  the 
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program    is    shown    in   Appendix   III    (labeled   Results   of   Data   Set 
II),    and   agrees  with   the   calculation   in  Appendix   IV. 


CHAPTER  V 


RECOMMENDATIONS 


At  the  onset  of  this  project,  an  objective  statement  was  set 
forth  which  was 

"To  develop  a  program  to  minimize  the  cost  of  natural 
gas  from  a  large  number  of  sources  over  a  long 
time  frame." 

The  results  presented  by  this  paper  satisfy  this  objective. 
However,  as  with  any  problem,  there  are  a  number  of  modifications 
which  could  be  made  to  make  the   program  more  realistic. 

The  greatest  room  for  improvement  in  the  solution  that  the 
author  sees  lies  with  the  problem  of  December  price  escalation. 
This  escalation  causes  flow  rates  to  be  assigned  to  some  sources 
in  an  order  that  is  not  chronological.  The  problem  with  this 
unordered  assignment  of  flow  rates  results  from  the  fact  that  the 
maximum  rate  allowed  any  source  is  decreased  after  each  flow  rate 
is  assigned.  As  a  result,  the  maximum  rate  and  reserves  for  some 
sources  are  not  in  the  proper  order.  As  an  example,  the  results 
from  Data  set  II  shown  in  Appendix  III  give  the  reserves  for 
source  one  as  569,520  in  time  step  9,  and  567,520  in  time  step  1. 
In  actuality,  the  reserves  must  decrease  with  time. 

Another  area  that  could  be  modified  is  the  method  used  to 
calculate  initial  reserves.  The  method  used  in  this  solution 
assumes  that  every  well  has  at  least  ten  years  reserves 
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available.  A  more  satisfactory  solution  would  result  if  actual 
well  data  could  be  analyzed  for  each  source  so  that  a  more 
realistic  number  could  be  used.  This  was  not  done,  primarily,  in 
an  effort   to   protect   proprietary   information   from   publication. 

The  method  used  to  decrease  maximum  allowable  flow  rates 
following  each  rate  assignment  could  also  be  improved.  The 
percentage  decrease  used  in  this  program  is  better  than  allowing 
the  maximum  rate  to  remain  constant,  but  the  author  believes  a 
better  method  exists.  This  is  an  area  that  could  be  the  subject 
of   further   investigation. 
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COMPUTER  PROGRAM  NUMBER  I 


DIMENSIO 
CHARACTE 
REAL  MIN 

mco  s»i  o 

INO*5 

OPR1*0. 

IONGPA=C 

COST-0. 

TMIN«0. 

LS*0 

NSCC=1 

1RECL*70 

ICCT»*1 

BESCAL*1 

DECESC=1 

0PEN(U, 

BESCAL-P 

OECESC-0 

MEXP-COU 

R(I)-DES 

INO-SOUP 

OPRI-PRI 

IONGPA-N 

NS-NUMBE 

NTS-NUMB 

MONTH- (A 

FACINF- 

ALL-ALLO 

MIN-MINI 

MAX-MAXI 

ISTCD-ST 

IFLDCO-F 

NAME-WEL 

2DPTM-IN 

PRKE1-P 

MCPACD- 

NGPACO 
NSU-NUMO 
COEFF-  C 
COST-COS 
RES-WELL 
TM1N-TOT 
LS-IDENT 
MONO-ORI 
NSCC-NUf 
NCS-NUMB 
ICCTR-CC 


N  R(60) 
R*2C  NAME 
/MAX 
0 


.005 
.0*5 

ACCESS*' 
ASE  CSCA 
ECEMBER 
NTER  INO 
IRED  RAT 
CE  FROM 
C?  OF  TH 
GPA  CODE 
R  CF  SOU 
ER  OF  TI 

NUMBER) 
IHFLATIO 
WAPLE 
MUM  FLOW 
MUM  FLOW 
ATE  CODE 
IfLD  COD 
L  NAM? 
DICATES 
RICE  OF 
f.GPA  COD 
«1   OTHE 
ER  OF  SO 
OEFFICIE 
T  OF  THE 

RESERVE 
AL  OF  MI 
I  FUNG  N 
GINAL  ST 
BER  OF  S 
EP  OF  SO 
EFFICIEN 


DIP' /RE CL=IRECL/RCDS=IRCDS/F OP M  =  " FORMATTED') 

LATION  FACTOR  (PER  MONTH) 

ESCALATION  FACTOR 

ICATING  THE  NUMBER  OF  DECEMBERS  THAT  HAVE  FASSED 

E  FOR  EACH  TIME  PERIOD 

WHICH  DATA  IS  READ 

E  HIGHEST  PRICED  GAS 

OF  THE  HIGHEST  PRICED  GAS 
RCES 
HE  STEPS 

THE  MONTH  OF  THE  FIRST  TIME  STEP 
N  FACTOR 

FOR  A  GIVEN  SOURCE 
FOR  A  GIVEN  SOURCE 


A  DEEP  (0)  OR  A  SHALLOW  (1)  WELL 

THE  GAS  IN  THE  SOURCE  BEING  CONSIDERED 

E  FOR  THE  SOURCE  BEING  CONSIDERED.   102  OP  108 

RWISE  NGPACD=0 

URCES  USED.   A  COUNTER 

NT  FOR  THE  COST  EQUATION  (  AN  ARRAY) 

GAS  OVER  ALL  TIME  PERIODS. (  TC  eE  MINIMIZED) 
S 

NIMUM  TA<ES 

UMP-ER  FOR  THE  HIGHEST  PRICED  SOURCE 
ARTIMG  MONTH 
ORTED  CARDS  COUNTED 
RTED  CARDS 
T  COUNTER 


READ  IN  NUMBER  OF  SOURCES/  NUMBER  OF  TIME  STEPS/  3EGINNING 
MONTH/  MONTHLY  INFLATION  FACTOR/  AND  DESIRED  MONTHLY  RATES 
AND  WRITE  TO  UNIT  17 

READCINO/1000)  NS/ NTS/ MONTH/ FACINF 
READCINO/2000)  (R(I)/I»1/NTS) 
»RITE (17,1 000) NS /NTS/MONTH/ FACINF 
WRITE(17/2G00)(P(I),I=1,NTS) 

SEARCH  FOR  HIGHEST  PRICED  GAS  AND  SAVE  INFO  ON  THAT  WELL/ 
DO  RESERVE  CALCULATIONS/  FIND  THE  SU*  OF  THE  MNIMUMS  AND 
WRITE  DATA  TO  A  RECORD 


00  80  I*1/NS 

READ (INC /JOCO)  ISTCD/IFLDCD/NAME/IDPTH/ALL/MIN/MAX/PRICE1/NGPACD 

IFUDPTH.SQ.1)  GO  TO  120 

RES»120.»ALL 

GO  TO  130 
120  RES»240.»ALL 
130  CONTINUE 

IF(PRICE1.LE.0PR1)  GO  TO  90 

0PR1»PRICE1 

L5=NS      SLA?EL  THIS  SOURCE  AS  THE  HIGHEST  PRICED  SOURCE 
90  TMIN^**NGI*C0    SS*VE  ™E  NGP*  C0I>E  °F  ™E  HIGHEST  BRiCED  SOURCE 
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COiir'UTiSR    PROGRAM    NUMBER    I 
CONTINUED 


c 

C      STORE  SCURCE  INFO  IN  A  DIRECT  ACCESS  FILE 
C 

WRITE  (U' 1/3000)  ISTCD/IFLDCD/KAME/IDPTH/RCS/MIN/MAX/PRICE1/NGPACD 
80  CONTINUE 
C 

C      WRITE  LS/0PR1/ICNGPA,  AND  TMIN  TO  POSITION  1  OF  THE  COEFFICIENT  FILE 
C 

WRITEC17,3500)LS/0PR1,I0NGPA/TriN 
C 

C      CALCULATION  OF  COEFFICIENTS 
C 

MCN0=MONTH 

NSU  =  1   ICOUNTER  FOR  NUMBER  OF  SOURCES  USED 
60  BEAD C14'NSU,310C)  PR  I C E 1 / NGP AC D 
M0NTH*M0N0 

IFINSU.F0.L5)  GO  TO  140   SELIMINATE  HIGHEST  PRICED  GAS 
DO  50  J»1,NTS 

IFCNGPACD.EO.DGO  TO  300  =!IS  THIS  SOURCE  NGPA  102  OR  103 
<FLAG=0 
GO  TO  310 
300  KFLAG*1 

310  IFdONGPA.eO.D  GO  TO  320  SIS  HIGHEST  PRICED  SCURCE  NGPA  1C2  OR  108 
LFcAG*0 
GO  TO  330 
320  LFLAG=1 
530  CONTINUE 

IFCKFLAG.NE.O.CP.LFLAG.NE.O)  GC  TO  340 

C0EFF=-CPR1*BESCAL*«CJ-1)/FACUF**(J-1)  +  PRICE1«SESCAL*«U-1)/ 
AFACINF««(J-1) 
GO  TO  460 
C 

C      DETERMINE  THE  NUMBER  OF  DECEMBERS  THAT  THE  STUDY  HAS  PASSED  THROUGH 
C 

340  KEXP*M0NTH/12 

IFCKFLAfi.NE.O)  GO  TO  410 

C0EFF«-0PR1»BE3CAL**(J-KEXP-1)»DECESC««ME'XP/FACINF««  (J-1)*PRICE1 
A»?ESCAL*«(J-1)/FACINF»*(J-1) 
GO  TO  460 
410  IF(LFLAG.EC.O)  GO  TO  420 

C0EFF=-CPR1*BESCAL*«(J-MEXP-1)*DECESC**MEXP/FACINF«*(J-1)«PPICE1 
A*9ESCAL»»(J-MEXP-1)«DECESC**(J-MEXP-1)/FACINF*«(J-1) 
GO  TO  460 
420  C0EFF  =  -CPR1*RESCAL.»*(J-1)/FACINF«*<J-1)*PRICE1«t)ESCAL**U-PEXP-1) 

A*DECESC»*.1EXP/FACINF«»(J-1) 
460  CONTINUE 

ICCTR*1CCTR+1 

TA<E«MIN  SINITIAL  RATE  FOR  EACH  IS  MIN 
M0NTH=M0NTH*1 
50  WRITt(15/4000)NSU/J/COEFF,TAKE 

140  IF(NSU.EQ.NS)  GO  TO  70 
NSU*NSU+1 

GO  TO  60 
70  CONTINUE 
C 

C      CALCULATE  INITIAL  COST  AND  WRITE  THAT  COST  TO  COEFFICIENT  FILE 
C 

MONTH*MONO 

DO  110  I«1/NTS 

IFCI0NGPA.I.E.1)  GO  TO  141 

MEXP»M0NTH/12 

C0ST  =  C0ST»0PR1»DECESC««MEXP*BESCA'L**<I-1-MEXP)/FACINF**<I-1)«RCI) 
GO  TO  110 

141  COST*C0ST*OPR1*PESCAL*»(I-1)/FACINF»*(I-1)*RCI) 
110  M0NTH*N0NTH*1 

C 

W»ITE(17,3700)  COST 
1000  FCRMATCI4,I2,I2,F5.3) 
2000  F0PMATC6E13.5) 

3000  F0RMAT(I2,I4,A20/I1,F10.0/F10.C/F10.0,F6.4,I1) 
3100  F0RMAT(57X/F6.4,I1) 
3500  F0R"!AT<I4,F6.4/I1,F15.0) 
3700  F0RMAT(F15.2) 
4G0U  F0RMATCI4,I2,F8.4,F10.0) 

STOP 

END 
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COMPUTER  PROGRAM  NUMBER  II 


DIMENSION    R( 

CHAIUCTER*20 

REAL  MIN/MAX 

IPC05-1 CO 

IN0»5 

OPH1»0. 

I0NGPA»0 

C0ST«0. 

TMIN-0. 

L3*0 

NSCC"1 

IRECL«7G 

ICCTR»1 

open<u,acce 
bexp-counter 
rci)-desiped 
ino-source  f 
opri-price  0 
i0n5pa-ngpa 
ns-number  of 
nts-number  0 
m0nth-(a  num 
min-minibub 
max-max;mum 
istcd-state 
ifldcd-field 
name-well  na 
idpth-indica 
pricc1-price 
ngpaco-  hgpa 

H5PACD«1 
NSU-NUMBER  0 
COST-COST  OF 
RES-WELL  RES 
TMIN-TOTAL  0 
LS-IDENTIFYI 
MONO-CRIGINA 
NSCC-f;UMBER 
NCS-NUMRSR  0 
ICCTR-COEFFI 


60) 
NAME 


SS='0IR' 
INDICAT 
RATE  FO 
ROM  WHIC 
F  THE  HI 
CODE  OF 

SOURCES 
F  TIME  S 
BER)  THE 
FLOW  FOR 
FLOW  FOR 
CODE 

CODE 
ME 

TES  A  DE 

OF  THE 

CODE  FO 

OTHERWIS 

F  SOURCE 

THE  GAS 

ERVES 

F  MINIMU 

NG  NUM9E 

L  STAPTI 

OF  SORTE 

F  SORTED 

CIENT  CO 


/RECL*IRECL,RCDS*IRCDS/FORM«' FORMATTED') 

ING  THE  NUP8ER  OF  DECEPBERS  THAT  HAVE  PASSED 

P  EACH  TIME  PERIOD 

H  DATA  IS  READ 

GHEST  PRICED  GAS 

THE  HIGHEST  PRICED  GAS 

TEPS 
MONTH  OF  THE  FIRST  TIME  STEP 
A  GIVEN  SOURCE 
A  GIVEN  SOURCE 


EP  (0)  OR  A  SHALLOW  (1)  WELL 
GAS  IN  THE  SOURCE  BEING  CONSIDERED 
R  THE  SOURCE  BEING  CONSIDERED.   102  OR  108 
E  NGPACD=0 
S  USED.   A  COUNTER 
OVER  ALL  TIME  PERIODS. (  TC  PE  MINIMIZED) 

M  TAKES 

R  FOR  THF  HIGHEST  PRICED  SCURCE 

NG  »ONTH 

D  CARDS  COUNTED 

CARDS 
UNTER 


READ  IK  NUMBER  OF  SOURCES,  NUMBER  OF  .TIME  STEPS/  BEGINNING 
MONTH,  MONTHLY  INFLATION  FACTOR/  AND  DESIPED  MONTHLY  RATES 


READC17/1000)  NS/NTS/MONTH/FACINF 
READ<17,2000)  <R(I)/I»1/NTS) 

NCS*(NS-1)«NTS 

READ  IN  HIGHEST  PRICED  SOURCE/  THAT  SOURCES  PRICE  AND  NG°A  CODE 
THE  SUM  OF  THE  BINS/  AND  THE  COST 

REA0(17/380C)LS/0PR1/I0NGPA,TMIN/C0ST 
WRITE(6/39C0)  COST 

DECREASE  THE  DESIRED  RATE  FOR  EACH  TIME  PERIOD  3Y  THE  SUM  CF 
THE  MINIMUM  TAKES 


DO  220  I»1/NTS 
R.(I)»»(I1-.TMIH   . 
IF(RCI)  .GT.O.O)  GO  TO  220 
WRITE(6/7000>  I 
GO  TO  9600 
220  CONTINUE 

BRING  IN  SORTED  COEFFICIENT  CARDS  ONE  AT  A  TIME 

160  1CCTR*ICCTR*1 

READ<16/OCC)NSU/J/COEFF/TAKE 

READ  SOURCE  CARD  CORRESPONDING  TO  SOURCE  NUMBER  NSU 

»EADCn*NSU/3300)ISTCD/IFLDCD/NAME,RE3/MIN/MAX 
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COMPUTER  PROGRAM  NUMBER  II 
CONTINUED 


c 

CALCULATE  THE  AMOUNT  OF  GAS  TO  TAKE  FROM  A  SOU6CE 

TAKE*«AX 

IF(RCJ)  .GT.0.0)  GO  TO  210 

TAKE*MIN 

GO  TO  211 

210  IF(TAifE.GT.PU))  TAKE'R  C  J  )  *MIN 
MAX*KAX-TAKE/RES«MAX 
RES«RES-TAKE 

211  CONTINUE 
C 

C      REDUCE  THE  DESIRED  RATE  BY  THE  AMOUNT  OF  THE  TAKE 

RU)«RCJ)-TAKE*MIN 
C 

C      PRINT  OUT  RESULTS 
C 

IFCNSCC.EQ.1)  WRITE<6,5000) 

WRITE<6,6000)  ISTCD,IFLDCD/NAME/NSU,J,COEFF,TAKE,RES 

C      WRITS  DATA  BACK  TO  SOURCE  CARD  DATA  BASE 

NCS*(NS-1)*NTS 

WRIT  ECU  'HSU/3300)  I  STCD/I  FLDCD/NAME,  RES/MIN/MAX 
C0ST«C0ST*COEFF«TA<E 
IF(NSCC.EO.NCS)  GO  TO  170 
NSCC*NSCC*1 
GO  TO  160 
170  CONTINUE 
C 

READ  INFO  ON  HIGHEST  PRICED  SOURCE  AND  ASSIGN  TAKES  TO  THAT  SOURCE 

READ(14'L3,3  3C0>  ISTCD, I FLD CD, NAME, RES /KIN," AX 
DO  190  I*1/NTS 
IFCR(I).LE.O.O)  GO  TO  200 
TAKE»«1AX 

IF(TAKE.GT.PO)  TAKE*R(I)+MIN 
MAX»MAX-TAKE/»ES»MAX 
RES«RES-TAKE 
R(I)«RCI)-TAKE+KIN 
IFCR(I).GT.O.O)  WRITC(6/9000)  I 

WRITE(U'I,3  300)  ISTCD,IFLDCO,KAME,RES,MIN/MAX 
WR1TEC6/950C)  ISTCD/ I FLDCD/NAHE/LS/ I/TAKE/RES 
GO  TO  190 
200  TAKE*MIN 

WPITEU/O5  30)  I STCD/IFLDCD/NAME/LS, I, TAKE/RES 
1'0  CONTINUE 

-RITECo/SOOO)  COST 
1000  F0R»ATU4,I2,I2,F5.3) 
2000  F0RNAT(6E13.5) 

3200  FORMATCI2,l4,»20,ix,M0.0,10X,F10.0> 
3  300  FORMAT<I2/U,A2C,1X,F10.0/F10.C,F10.0) 
3800  F0RNATCI4,F6.4,I1,F15.0,/F15.2) 
3900  FORMATCF15.2) 

4000  FORMAT(I4,I2,F8.4,F10.0) 
5000  " 


6000 
7000 


FORMATdX,  'STATC1X,'  FIELD  ',1CX,'  NAME  ' / 1  OX/ • NSU ' / 2 X ,  -PPR  ICC  • 
A4X/.C0EFF./5X, -TAKE' /10X/ -PES', /,1X/'C0DE'/2X/ -CODE') 
.«?rnI(2X'i:'3X'I4'2X'A20'2X'U'3X'I2'4X'F8'4'2'</F10.0/2X/ 

FORMAT<//, 'DESIRED  RATE  FOR  TIME  PERIOD  ',11, •  WAS  LESS  THAN'/ 
A'THC  MINIMUM  REQUIRED  FLOW  RATE')  * 

3000  F0RMAT('-',2X,'THE  FINAL  COST  CF  THE  GAS  IS  '/F15.2) 
9<J00  FORMATd/,'  DESIRED  RATE  FOR  TIME  PERIOD  ',12,'  COULD  NOT  SE 

nil  c0o^;NTur':2'3x'u'2x'A20'2x'u'3x'i2'ux'n°:c'"'"o-°' 


»ET') 


STOP 

DEBUG  INIT/UNITC21 ) 

END 
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FLOW  CHART  FOR  PROGRAM 
NUMBER  I 
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Read:Number  Of  Sources, 

Number  Of  Time  Steps, 

Beginning  Month, 

Inflation  Factor, 

Monthly  Desired  Rates 


Write  Data  Above  To  File 

17  For  Reference  By 
Program  II 


NO 


■> 


Read  Data  For  One  Source 
At  A  Time 


X 


Calculate  Reserves  Based 
On  Depth 


± 


YES 


Does  This  Source  Have  The  Highest]   j  Save  Source 
Price  Encountered  Thus  Far?   J— ^Number,  Price 

And  NGPA  Code 


NO 


<r 


Increase  Total-Minimum 
By  This  Minimum 


Write  This  Sources 
To  Source  F 


Information  1 
ile | 


1 


-fis  This  The  Last  Source?] 

£YES 


Write,  For  Source 

With 

Highest  Cost: 

Source 

Number 

,  Price,  and 

NGPfl 

Code 

Flag 

and 

Total  Of 

Minimum  Fl 

.ow  Rates 

To  File 

17  For 

Reference  By 

Proqram  II 

\ 


CONTINUED 


FLOW  CHART  FOR  PROGRAM  NUMBER  I 
CONTINUED 


39 


V 


Set  Variable  Mono  To  The 
Beginning  Month 


Set  ICCTR  And  NSU  Counters  To  1 


I 


I  For  Source  #NSU,  Read: 
*yj   Price  And  NGPA  Code 


J 


jSet  MONTH  Equal  To  MONO 


Is  This  The  Highest  Priced  Source?- 


NO 


^l 


Calculate  Coefficients  Considering 
NGPA  Code  For  This  Source  And 
For  The  Highest  Priced  Source_ 


NO 


t 

Increment 

ICCTR 

By 

1 

i 


Set  Take  Initially  To 
Minimum  Allowed 


1 


Write :Number  Of  Source, 

Time  Period  Being  Considered, 

Calculated  Coefficient, 

And  Take  To  Coefficient 

File 


Is  This  The  Last  Time  Step? 
!   YES 


YES 


Incrementj^- 
NSU 


NO 


Is  This  The  Last  Source?) 


V 


YES 
CONTINUED 
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FLOW  CHART  FOR  PROGRAM  NUMBER  1 
CONTINUED 


Jl 


Ca 

lculate  The 

Initial 

Cost  i 

Based  On  The  Highe 

»t 

Priced 

Source 

| 

I 


Write 

Initial 

Cost 

To  File 

17 

For  Ref 

=rence 

By 

Prog 

ram  II 

■ 

FLOW  CHART  FOR  PROGRAM 
NUMBER  II 
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Read  From  File  17:  Number  Of 
Sources,  Number  Of  Time 
Steps,  Starting  Month, 
Inflation  Factor,  And 
Desired  Rates 

v 

Set  The  Number  Of  Coefficients 
To  Be  Sorted  at 

(NS-1) *NTS 


J 


Read  In  Highest  Priced  Source 
Information  From  File  17 


Decrease  The  Desired  Rate  For 

Each  Time  Step  By  The  Sum 

Of  The  Minimums 


YES 


Are  Any  Of  The  New  Desired 
Rates  Less  Than  Zero 


y  NO 

Read  Coefficient  Cards  In 

The  Order  That  They 
Appear  In  The  Sorted 
Coefficient  File 


>J 


Print  Out 
Error  Message 
And  End 


Read  Source  Card  Corresponding] 
To  Coefficient  Card  Just  Read| 


Calculate  The  Maximum  Takej 
Considering  The  Reserves,* 
Maximum  Allowable  Rate, 
And  Desired  Rates 


CONTINUED 
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A 


FLOW  CHART  FOR  PROGRAM  NUMBER  II 
CONTINUED 


1 


Reduce  The  Desired  Rate  For  The  Time 

Step  Represented  By  The  Last 

Coefficient  By  An  Amount  Equal 

To  The  Take  Assigned 


JL 


IPrint  Out  The  Results  For  This 
Source/Time  Step  Combination 


Write  Information  Back  To  The  Source  Card 
File,  Including  Changed  Values  For 
The  Source  Reserves  and  Maximum 

Allowable  Flow  Rate       


Adjust  The  Total  Cost  Term 
For  This  Take 


Increment  # 
Of  Sorted 
Coefficients 
Counted 


NO 


Have  All  Sorted  Coefficients 
Been  Considered? 


J_ YE 


YES 


Read  Information  On  The  Highest 
Priced  Source,  And  Assign 
A  Take  For  Each  Time  Step 


Print  Out  Each  Source/Time  Step 
Results 


I 


Print  Out  Final  Cost 


APPENDIX    III 
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DATA    SET    I 


44 


CC7UG091.01C 
4.98060E  06 
4.93G60E  06 
4.98060E  06 
4.93060E  06 
4.9J060E  06 
4.98060E  06 
4.9306CE  06 
06 


4.98060E  C6 

4.98060E  06 

4.98060E  06 

4.9S060E  06 

4.98060E  06 

4.95060E  06 

4.9S060E  06 

4.98060E  06 

4.98060E  06 

4.95060E  06 


4.98060E 

4.950oOE  06 

4.9S06QE  06 
0110C0FIJL0  ON: 
021100FISL0  TWO 
0312COFIELO  THREE 
031300FIELO  FOUR 
031400FIEL0  FIVE 
G21500FIELO  SIX 
0316G0FIELD  SEVEN 
O11700FIEL0  EIGHT 
011300FIELO  NINE 
021900FIEL0  TEN 
O32OOGFIEL0  ELEVEN 
022100FIEL0  TWELVE 
0422C0FIEL0  THIRTEEN 
032300FIELD  FOURTEEN 
0324C0FIELD  FIFTEEN 
022500FIELD  SIXTEEN 
032600FIEL0  SEVENTEEN 
02270CFIELD  EIGHTEEN 
0323CGFIELD  NINETEEN 
C22900FIEL0  TWENTY 
033000FIEL0  TWENTY-ONE 
0231C0FIEL0  TWENTY— TWO 
013200FIELO  TWENTY-THREE 
013300FIELO  TWENTY-FOUR 
043400FIEL0  TWENTY-FIVE 
053500FIELO  TWENTY-SIX 
013600FIELO  TWENTY-SEVEN 

0237COFIELO  twenty-eight 

023800FIELO  TWENTY-NINE 
023900FIEL0  THIRTY 
034000FIEL0  THIRTY-ONE 
0241C0FIELD  THIRTY-TWO 
034200FIELO  THIRTY-THREE 
024300FIELD  THIRTY-FOUR 
034400FIEL0  THIRTY-FIVE 
044500FIEL0  THIRTY-SIX 
064600FIEL0  THIRTY-SEVEN 
044700FIEL0  THIRTY-EIGHT 
0243C0FIEL0  THIRTY-NINE 
014900FIELD  FORTY 
025000FIELO  FOPTY-ONE 
015100FIEL0  FORTY-TWO 
025200FIEL0  FORTY-THREE 
3453C0FIELD  FCRTY-FOUR 
025400FIELO  FORTY-FIVE 
0155O0FIELD  FORTY-SIX 
025600FIELO  FORTY-SEVEN 
025700FIELD  FORTY-EIGHT 
0153CCFIELO  FCRTY-NINE 

O159C0FIELO  FIFTY 
066000FIEL0  FIFTY-ONE 
0461C0FIEL0  FIFTY-TWO 
0162C0FIEL0  FIFTY-THREE 
046300FIEL0  FIFTY-FOUR 
01640QFIEL0  FIFTY-FIVE 
0265C0FIELD  FIFTY-SIX 
0166C0FIEL0  FIFTY-SEVEN 
0667C0F:eLC  FIFTY-EIGHT 
026300FIELD  FIFTY-NINE 
Q2690CFIcLD  SIXTY 
027OC^FIELD  SIXTY-ONE 
G571CuF1EU0  SIXTY-TWO 
0672C0'IEL0  SUTY-THPEE 

0273c:Fiei.r  sixty-four 

G27iC0FIELCi  3KTY-=IVE 
O27500f;el0  SIxty-six 
^47^,C:FIELD  SIXTY-SEVjN 
077700FIEL0  SIXTY-EIGHT 
0673C;eIEL0  SIXTY-NINE 
0o7?:CFIcL0  SEVENTY 

os^oooFiei:  seventy-cn= 


98060E  06 

98060E  06 

"3060E  06 

9306OE  06 

"3060E  06 

98060E  06 

95060E  06 

93060E  06 

4.98060E  06 

4.93060E  06 

2373. 

2888. 

1769<>7. 

1605205. 

300653. 

8238. 

13430. 

30600. 

25S19. 

521858. 

25819. 

46967. 

2038. 

155000. 

T1CO0O. 

113893. 

44334. 

162134. 

145403. 

7C493. 

»°1S. 

66077. 

2543. 

4841. 

110"21 . 

769649. 

1359. 

178-4. 

104041. 

18770. 

202222. 

7644. 

3058. 

49685. 

10786. 

2463. 

8758. 

42«75. 

8153. 

16477. 

6E200. 

36775. 

340. 

034. 

54951. 

255. 

11466. 

25564. 

34822. 

2373. 

1529. 

50789. 

2633. 

10192. 

2203. 

6235. 

<07. 

32359. 

13079. 

107=o. 

34567. 

IIS1. 

6370 , 

101O;. 

15283. 

32650. 

10701  . 

64*o;i. 

30575. 

1 •40000. 

1 4  »  6 0  . 


92060 
98060 
98060 
98060 
95060 
98060 
98060 
93060 
98060 
98060 

1000. 
300. 

10C0. 

1000. 

10CO. 

1000. 

1000. 

1000. 

1000. 

1000. 
SCO. 
30C. 

1000. 

1000. 

1000. 

1000. 

1000. 

10CC. 

10CO. 
6C0. 
60C. 

1000. 

1000. 

1000. 

1000. 

1000. 

1000. 
300. 

1000. 
4C0. 

1000. 

10CC. 

1000. 
400. 
100. 

1000. 

1000. 

10C0. 

1000. 
500. 

1000. 

10CO. 
100. 
917. 
900. 
1C0. 

10C0. 
30C. 
7C0. 

1000. 
1000. 

100:. 

10C0. 
1000. 
1000. 

3C0. 

100. 
1000. 
1000. 

400. 

100. 

1000. 
1000. 

1000. 
1000. 
10CC. 

loec. 

1000. 
'000. 
1000. 
1000. 


06 
06 
06 
06 
06 
06 
06 
06 
06 
06 


4. 

4. 

4. 

4. 

4. 

4. 

4. 

4. 

4. 

4. 

2567 

2833 

173666 

1575000 

295C00 

3053 

18083 

2817123 

514000 

25333 

33C00 

5O000 

2200 

152083 

304167 

122925 

43500 

216511 

232500 

117000 

8750 

594685 

5C00 

4750 

108833 

195Cr00 

1344 

4247 

114750 

24000 

200000 

7500 

3000 

244665 

10583 

2417 

9253 

46383 

10000 

2C750 

67032 

36C33 

333 

917 

155C00 

275 

11250 

461221 

163023 


06 
06 


06 
06 


O8060E  06 

9806CE  06 

93060E 

9506CE 

9e060E  06 

98060E  06 

9806CE  06 

93060E  06 

98060E 

98060E 

.43340 

.47620 

.49870 

.57800 

.61361 

.65690 

.67210 

.74510 

.7?990 

.74691 

.95310 

1.10O40 

1.12450 

1.16461 

1 .38900 

1.49500 

1.61960 

1.67650 

1.77440 

.1.83110 

.2.10550 

.2.15861 

.2.22831 

.2.33000 

.2.35930 

.2.40000 

.2.45370 

.2.6C330 

.2.67510 

.2.74750 

.2.75300 

.2.80001 

.2.86830 

•2.9C200 

.2.94710 

.2.00230 

.3.03980 

.3.08660 

.3.12000 

.3.16431 

.3.20950 

.3.21141 

.3.22910 

.3.23290 

.3.29840 

.3,30001 

.3.35601 

.3,41650 

.3.46030 


3C00 

1650 
54817 

3000 
11C00 

3C00 
11000. 

1C0O. 
34925. 

13000. 

104  16. 

35500. 

1467. 

6875. 

1 1000. 

lonoo. 

34200. 

1H50. 

45»0°0. 

3CC00. 

4°6?P00. 

14553. 


3.53711 
3.53711 
3.5946-1 
.3.65001 
.3.69250 
.3.73221 
.3.76501 
.3.80821 
.3.53931 
.3.96001 
.4.04171 
.4.05001 
.4.07601 
.4.15751 
4.270C1 
4  . 4srni 

4.52001 
4.573X1 

4. TOCO 
4.»162n 
5.1C11" 
•».  nnonn, 


4.9S060E  06 

4.98060E  C6 

4.98060E  06 

4.98060E  C6 

4.9S06CE  06 

4.98060E  06 

4.9S06CE  06 

4.9506CE  06 

4.9S060E  C6 

4.98060E  C6 
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DESIRED  RATE  SUM  FOR  ALL  FIELDS  IN  TIME  PERIOD  ONE 

DATA  SET  I 

Field  Assigned  Rate       Field  Assigned  Rate 

2567.             37  1000. 

38  1000. 

39  1000. 

40  500. 

41  1000. 

42  1000. 

43  100. 

44  917. 

45  900. 

46  100. 

47  1000. 

48  300. 

49  700. 

50  1000. 

51  1000. 

52  1000. 

53  1000. 

54  1000. 

55  1000. 

56  300. 

57  100. 

58  1000. 

59  1000. 

60  400. 

61  100. 

62  1000. 

63  1000. 

64  1000. 

65  1000. 

66  1000. 

67  1000. 

68  1000. 

69  1000. 

70  1000. 

71  1000. 
36                 1000.          Total  Rate  4980600. 


2 

2833 

3 

173666 

4 

1575000 

5 

295000 

6 

8083 

7 

18083 

8 

2817128 

9 

12490 

10 

25333 

11 

800 

12 

800 

13 

1000 

14 

1000 

15 

1000 

16 

1000 

17 

1000 

18 

1000 

19 

1000 

20 

600 

21 

600 

22 

1000 

23 

1000 

24 

1000 

25 

1000 

26 

1000 

27 

1000 

28 

300 

29 

1000 

30 

400 

31 

1000 

32 

1000 

33 

1000 

34 

400, 

35 

100 
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DATA    SET    II 


0C1220041.O10 

2.CCO00E    03      1.90C0CE    03  1.*0000E    03 

1.40000c    03       1.30000:    03  1.200008    03 

1.10000E    03      1.100005    03  1.10000E    03 

.     1.1000CE    03      1.10000=    03 

011000FIELO  ONE           1  2378, 

03HC0FIELD  FIVE           1  300653. 

032000FIELD  ELEVEN         0  25819. 

02250CFIELD  SIXTEEN        0  113?93. 

023100FIELC  TWENTY-TWO     0  66077. 

013600FIEL3  T'wE!4TY-3EVEN   0  135'. 

O3-.2O0FIELO  THIRTY-THREE   0  3058. 

0447Q0FIEL0  THRITY-EIGHT   0  42975. 

0453C0FIEL0  FCRTY-FCUR     0  934. 

0153C0FIEL0  FORTY-NINE     0  34«22. 

016500FIEL0  FIFTY-SIX      0  2203. 

02710QFIEL0    SIXTY-TWO               0  345t7. 


1.7C000E  03 
1.1C090I  03 

1 .1CC00E  03 


1.60O0C5  03 
1.1000CS  C3 
1.100005  03 


1.50000!  03 
1 . 1C0CCE  C3 
1.10000E  03 


1030. 

2567. 

,43  3  40 

ccoc. 

000000. 

.61361 

ccc. 

00000. 

.95? 10 

cooc. 

C0CC000.1 

.-«500 

oooc. 

000000.2 

.15361 

ooco. 

0000.2 

.45370 

0000. 

0000.2, 

.  ?6*30 

oooc. 

OOOCO .  3 

.0E660 

000. 

000.3 

.2*2*0 

coo. 

0C00CO.3 

.46030 

0000. 

0000.3 

.73221 

coco. 

00000.4 

.05001 

RESULTS  OF  DATA  SET  II 
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STATE  FIELD 

NA* 

CODE 

CODE 

1 

10C0 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

1 

1000 

FIELD 

ONE 

3 

1400 

FIELD 

FIVE 

3 

1400 

FIELD 

FIVE 

3 

1400 

FIELD 

FIVE 

3 

1400 

FIELD 

FIVE 

3 

1400 

FIELD 

FIVE 

3 

1400 

FIELD 

FIVE 

3 

1400 

FIELD 

FIVE 

3 

1400 

FIELD 

FIVE 

3 

1400 

FIELD 

FIVE 

3 

1400 

FIELD 

FIVE 

3 

2000 

FIELD 

ELEVEN 

3 

2000 

FIELD 

ELEVEN 

3 

1400 

FIELD 

FIVE 

3 

2000 

FIELD 

ELEVEN 

3 

2000 

FIELD 

ELEVEN 

3 

1400 

FIELD 

FIVE 

3 

2000 

FIELD 

ELEVEN 

3 

2000 

FIELD 

ELEVEN 

3 

2000 

FIELD 

ELEVEN 

3 

2000 

FIELD 

ELEVEN 

3 

2000 

FIELD 

ELEVEN 

3 

2000 

FIELD 

ELEVEN 

3 

2000 

FIELD 

ELEVEN 

3 

1400 

FIELD 

FIVE 

3 

2000 

FIELD 

ELEVEN 

3 

2OC0 

FIELD 

ELEVEN 

3 

2000 

FIELD 

ELEVEN 

3 

2000 

FIELD 

ELEVEN 

3 

2000 

FIELD 

ELEVEN 

3 

2000 

FIELD 

ELEVEN 

3 

2000 

FIELD 

ELEVEN 

3 

14  00 

FIELD 

FIVE 

3 

2000 

FIELD 

ELEVE;.- 

7 

2 

2000 

FIELD 

ELEVEN 

3 

noo 

FIELD 

FIVE 

3 

1400 

FIELD 

FIVE 

3 

14C0 

FIELD 

FIVE 

3 

1400 

FIELD 

FIVE 

NSU   PERIOD 

COEFF 

TAKE 

1     9 

-3.6307 

1200. 

1     1 

-3.6166 

200U. 

1    10 

-3.6128 

1100. 

1     2 

-3.5987 

1900. 

1    11 

-3.5949 

1100. 

1     3 

-3.5809 

1800. 

1    12 

-3.5771 

1100. 

1     4 

-3.5632 

1700. 

1    13 

-3.5594 

1100. 

1     5 

-3.5455 

160C. 

1    14 

-3.5418 

1100. 

1     6 

-3.5280 

1500. 

1    15 

-3.5242 

1100. 

1     7 

-3.5105 

1400. 

1    16 

-3.5068 

1100. 

,1     8 

-3.4931 

1300. 

1    17 

-3.4894 

1100. 

1    18 

-3.4721 

1100. 

1    19 

-3.4550 

1100. 

1    20 

-3.4379 

1100. 

2     1 

-3.4364 

0. 

2     2 

-3.3919 

0. 

2     3 

-3.3466 

0. 

2    4 

-3.3003 

0. 

2     5 

-3.2530 

0. 

2     9 

-3.2487 

0. 

2     6 

-3.2048 

0. 

2    10 

-3.1969 

0. 

2     7 

-3.1556 

0. 

2    11 

-3.1439 

0. 

3     9 

-3.1265 

0. 

3    10 

-3.1110 

0. 

2     8 

-3.1C52 

0. 

3    11 

-3.0956 

0. 

3     1 

-3.0919 

0. 

2    12 

-3.0897 

0. 

3    12 

-3.0803 

0. 

3     2 

-3.0766 

0. 

3    13 

-3.0650 

0. 

3     3 

-3.0614 

0. 

3    14 

-3.0498 

0. 

3     4 

-3.0462 

0. 

3    15 

-3.0347 

0. 

2    13 

-3.0342 

0. 

3     5 

-3.0311 

0. 

3    16 

-J. 0197 

0. 

3     6 

-3.0161 

0. 

3    17 

-3.0048 

0. 

3     7 

-3.0012 

0. 

3   ie 

-2.9890 

0. 

2     ? 

-2.9863 

0. 

2    14 

-2.9774 

0. 

3     19 

-2.9751 

0. 

3    20 

-2.3604 

0. 

2    15 

-2.9192 

0. 

2    16 

-2.3595 

c. 

2    17 

-2.7<??3 

0. 

2    15 

-2.7356 

0. 

RES 

569520. 
567520. 
56642C. 
564520. 
563420. 
561620. 
560520. 
55382C. 
557720. 
556120. 
555020. 
553520. 
552420. 
551020. 
54O"20. 
549620. 
547520. 
546420. 
545320. 
544220. 
72157920. 
72157«20. 
72157920. 
72157920. 
72157920. 
72157920. 
72157920. 
72157920. 
72157920. 
7215792C. 
3C98280. 
3C98280. 
72157920. 
3093280. 
3C98280. 
7215792C. 
3C932E0. 
3098280. 
3C9828C. 
3C982E0. 
3C982e0. 
3C9B280. 
3098280. 
72157320. 
3098250. 
3C9S2E0. 
309S2E0. 

3C98260. 

3C9S2S0. 

3098260. 

JC°82?d. 
721579Z0. 

JC98280. 

3C9o280. 
72157920. 
72157920. 
72157O20. 
7215792C. 
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RESULTS  OF  DATA  SET  III 


FIELD  T*:rTY-Tn-?Zr 

FIELD  TriI?TY-T-isEE 

FIELD  THIRTY-THREE 

FIELD  THIS  TV -THREE 

FIELD  THIRTY-THREE 

FIELD  THIRTY-THREE 

FIELD  THIRTY-THREE 

FIELD  THIRTY-THRE~E 

FIELD  THIRTY-THREE 

FIELD  ThENTY-TWO 

FIELD  ThIRTY-THREE 

FIELD  THIRTY-THREE 

FIELD  THIRTY-THREE 

FIELD  TWENTY-TWO 

FIELD  THIRTY- THREE 

FIELD  THIRTY-THREE 

FIELC  THIRTY-THREE 

FIELD  THIRTY-THREE 

FIELD  THIRTY-THREE 

FIELD  THIRTY-THREE 

FIELD  THIRTY-THREE 

FIELD  THIRTY-THREE 

FIELD  TWEiJTY-TWO 


7 

9 

-1  .2a06 

1200. 

365760. 

7 

10 

-1  .2?42 

1 1  CO. 

364660. 

7 

11 

-1.2779 

110C. 

363560. 

f 

12 

-1.2715 

11CC. 

362460. 

7 

*  T 

-1  .2653 

1100. 

36136C. 

7 

1  u 

-1 .2590 

11CC. 

360260. 

7 

1  c 

-1  .252' 

110C. 

359160. 

7 

-I  < 

-1.2466 

1100. 

358060. 

7 

17 

-1.24.04 

11CC. 

356960. 

c 

3 

-1.2331 

c. 

7929240. 

7 

1  : 

-1.2342 

1100. 

355360. 

7 

1  9 

-1  .2281 

1100. 

354760. 

7 

2  31 

-1 .2221 

11CC. 

353660. 

- 

7 

-1.2026 

0. 

7929240. 

7 

1 

-1  .1317 

2000. 

351660. 

7 

2 

-1 .1759 

1900. 

349760. 

7 

T 

-1 .1700 

•  1300. 

347960. 

7 

i 

-1 .1642 

1700. 

346260. 

7 

5 

-1.1525 

1600. 

344660. 

7 

£ 

-1.1527 

1500. 

343160. 

9 

7 

-1 .1470 

1400. 

341760. 

7 

? 

-1  .1414 

1300. 

340460. 

; 

1C 

-1.1062 

C. 

7929240. 

RESULTS  OF  DATA  SET  IV 


FIELD 

SIXTY-TWO 

12 

1 

FIELD 

SIXTY-TWO 

12 

2 

FIELD 

SIXTY-TWO 

12 

3 

FIELD 

SIXTY-TwO 

12 

u 

FIELD 

SIXTY-TWC 

12 

5 

FIELD 

SIXTY-TWO 

1? 

6 

FIELD 

SIXTY-TWO 

12 

7 

FIELD 

SIXTY-TWO 

1  2 

E 

FIELD 

SIXTY-TfcO 

12 

9 

FIELD 

SIXTY-TWO 

12 

10 

FIELD 

SIXTY-TWO 

17 

11 

FIELD 

SIXTY-TWO 

12 

12 

FIELD 

SIXTY-TWO 

17 

13 

FIELD 

SIXTY-TWO 

1? 

14 

FIELD 

SIXTY-TWO 

12 

15 

FIELD 

SIXTY-TWO 

12 

16 

FIELD 

SIXTY-T».0 

12 

17 

FIELD 

SIXTY-TWO 

12 

18 

FIELD 

SIXTY-T*0 

12 

19 

FIELD 

SIXTY-TO 

1  2 

20 

2000. 
1900. 
130C. 
1700. 
1600. 
1500. 
1400. 
1300. 
1200. 
1100. 
1100. 
1100. 
1100. 
1100. 
1100. 
1100. 
1100. 
1100. 
1100. 
1100. 


4146040. 
4144140. 
4142340. 
4140640. 
4139040. 
4137540. 
4136140. 
4134340. 
4133640. 
4132540. 
4131440. 
4130340. 
4129240. 
412314C. 
4127040. 
4125940. 
4124340. 
4123740. 
4122640. 
4121540. 


FINAL  COSTS  FOR  DATA  SETS  II,  III,  &  IV 


Data  Set  II 

THE  FINAL  COST  OF  THI  3*3  I?  1 10Z3.5H 


Data  Set  III 

THE  PIHAL  CCST  Cr  TnE  «AS  IS         72*39.6 


Data  Set  IV 

THE  FINAL  CC3T  OF  THE  043  IS        105C23.67 


50 


APPENDIX    IV 
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MAXIMUM  FLOW  RATE  REDUCTION 
CALCULATIONS 


USING  THE  LARGEST  DATA  SET  (60*71), 
FIELD  1    TIME  PERIOD  1 

MAXIMUM  FLOW  RATE  ALLOWABLE  =  256  7.0 
FLOW  RATE  ASSIGNED  BY  PROGRAM  =  2567.0 

INITIAL  RESERVE  CALCULATIONS: 

240*ALLOWABLE=240*23  78=570720 

%  DECREASE  IN  RESERVES  CALCULATION: 

256 7/570720=4. 4978*10"3 

%  DECREASE  IN  MAXIMUM  ALLOWABLE  FLOW  RATE=4 . 4978*10"3 

NEW  MAXIMUM  ALLOWABLE  CALCULATIONS: 

MAX=MAX(1-4.4978*10~3) 

=  2567.0  (1-4.4978*"3) 

=2555.45 

THIS  RESULT  CHECKS  WITH  THE  MAXIMUM  ALLOWABLE  FLOW  RATE 
LISTED  BY  THE  PROGRAM  FOR  SOURCE  1  IN  TIME  PERIOD 
2. 


53 


COEFFICIENT  CALCULATIONS 


THIS  SAMPLE  CALCULATION  USES 
AND  CALCULATES  THE  COEFFICIENT 
ELEVENTH  TIME  PERIOD. 


THE  SMALLEST  DATA  SET  (21*12), 
FOR  THE  FIRST  FIELD  IN  THE 


FIELD  ONE  HAS  AN  NGPA  CODE  THAT  IS  NOT 
THEREFORE,  THE  PRICE  OF  THE  GAS  FROM  THAT  SOURCE  IS 
AN  EXTRA  FOUR  PERCENT  IN  DECEMBER. 


102  OR  108. 
NOT  ESCALATED 


FIELD  SIXTY-TWO  HAS  AN  NGPA  CODE  THAT  IS  102  OR  108. 
THEREFORE,  THE  PRICE  OF  THE  GAS  FROM  THAT  SOURCE  IS  ESCALATED  AN 
EXTRA  FOUR  PERCENT  IN  DECEMBER. 

SINCE  THIS  PARTICULAR  STUDY  BEGAN  IN  MONTH  FOUR,  OR  APRIL, 
MONTH  ELEVEN  REPRESENTS  FEBRUARY,  AND  THE  STUDY  HAS  PASSED 
THROUGH  A  DECEMBER. 


COEFF=     PRESENT  WORTH  OF  THE 
-  COST  ASSOCIATED  WITH 
THE  HIGHEST  PRICED 
WELL 


PRESENT  WORTH  OF  THE 
COST  ASSOCIATED  WITH 
THIS  WELL 


i — (initial  cost  of  highest  priced  source) 

(——(escalation  factor  for  months  other  than  Dec) 
I — (escalation  factor  for  December  past) 
i — (initial  cost  of  source  1) 

(escalation  factor  for  all 
months) 


c 


1 


COEFF=      -(4.0500) 


^  + 

(1.005)9    *     (1.045)1    +    .4334     *     (1.005) 


10 


(1.01) 


10 


(1.01) 

t 


10 


-(inflation   factor   for   all   months) 


C0EFF=  -3.5949 
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In  the  natural  gas  industry,  profitable  production 
depends  upon  many  factors.  These  factors  can  include  such 
things  as  the  cost  of  the  gas  from  each  source,  the  rate  of 
increase  of  that  cost,  the  applicability  of  take-or-pay,  and 
the  desire  for  rateable  production.  This  study  considered  one 
of  these,  namely  the  problem  of  cost  minimization.  In 
addition  to  minimizing  cost,  this  study  was  structured  such 
that  the  resulting  computer  program  was  capable  of  handling  a 
large  number  of  sources  through  many  time  steps.  These 
additional  problem  reguirements  made  this  study  different  from 
other  small  scale  cost  minimization  problems. 

Initially,  a  thorough  literature  search  was  completed. 
The  search  showed  that  very  few  problems  of  this  type  had  been 
solved.  Those  sources  produced  by  the  literature  search, 
dealt  with  problems  that  were  magnitudes  smaller  than  the 
problem  addressed  by  this  paper.  Without  exception,  the 
articles  found  presented  a  matrix  solution.  Because  the 
solution  to  this  problem  was  to  be  capable  of  handling  up  to 
5000  sources,  and  up  to  60  time  steps,  the  use  of  matrices  in 
the  solution  was  unacceptable.  A  solution  requiring  less 
active  computer  space  as  well  as  less  execution  time  was 
required. 

The  solution  obtained  required  sequential  use  of  two 
programs.  The  first  program  developed  a  coefficient  file  that 
contained  one  element  for  every  source/time  step  combination. 
This  file  could,  therefore,  contain  up  to  300,000  unknowns. 


The  coefficient  file  was  sorted  using  a  sort  processor  after 
the  first  program  had  completed  execution.  The  second  program 
processed  data  from  the  sorted  coefficient  file  and  determined 
appropriate    flow   rates. 

During  the  execution  of  the  second  program,  such  things 
as  variable  maximum  flow  rates,  well  reserves,  and  desired 
flow  rates  for  a  given  time  step  were  considered.  Each  of 
these  variables  was  constantly  changing  throughout  the 
execution   of    the   program. 

The  problem  for  which  this  program  was  developed  was 
except  in  magnitude,  a  typical  linear  programming  problem. 
However,  because  of  the  magnitude  of  the  problem,  the 
solution,    in   form,   was   atypical. 
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